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The Asymmetric Simple Exclusion Process (ASEP) plays the role of a paradigm in Non- 
Equilibrium Statistical Mechanics. We review exact results for the ASEP obtained by Bethe Ansatz 
and put emphasis on the algebraic properties of this model. The Bethe equations for the eigenvalues 
of the Markov Matrix of the ASEP are derived from the algebraic Bethe Ansatz. Using these equa- 
tions we explain how to calculate the spectral gap of the model and how global spectral properties 
such as the existence of multiplets can be predicted. An extension of the Bethe Ansatz leads to an 
analytic expression for the large deviation function of the current in the ASEP that satisfies the 
Gallavotti-Cohen relation. Finally, we describe some variants of the ASEP that are also solvable by 
Bethe Ansatz. 
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I. INTRODUCTION 

Equilibrium statistical mechanics tells us that the probability of a microstate of a system in equilibrium with a 
thermal reservoir is given by the Boltzmann-Gibbs law : if 7i is the Hamiltonian of the system, the probability 
distribution over the configuration space is proportional to exp(— /37i) where (3 is the inverse of temperature. This 
canonical prescription is the starting point for any study of a system in thermodynamic equilibrium : it has provided 
a firm microscopic basis for the laws of classical thermodynamics, has allowed us to describe various states of matter 
(from liquid crystals to superfluids) , and has led to a deep understanding of phase transitions that culminated in the 
renormalisation group theory. 

For a system out of equilibrium, the probability of a given microstate evolves constantly with time. In the long 
time limit such a system may reach a stationary state in which the probability measure over the configuration space 
converges to a well-defined and constant distribution. If the system carries macroscopic stationary currents that 
represent interactions and exchanges of matter or energy with the external world, this stationary distribution is 
generically not given by the canonical Boltzmann-Gibbs law. At present, there exists no theory that can predict the 
stationary state of a system far from equilibrium from a knowledge of the microscopic interactions of the elementary 
constituents of the system amongst themselves and with the external environment, and the dynamical rules that 
govern its evolution. The search for general features and laws of non- equilibrium statistical mechanics is indeed a 
central topic of current research. 

For systems close to thermodynamic equilibrium, linear-response theory yields the fluctuation-dissipation relations 
and the Onsager reciprocity relations. However, these properties of non-equilibrium statistical mechanics are only 
valid in the vicinity of equilibrium. Another route to explore and discover the characteristics of a general theory for 
systems out of equilibrium is through the study of mathematical models. These models should be simple enough to 
allow an exact and thorough mathematical analysis, but at the same time they must exhibit a complex phenomenology 
and must be versatile enough to be of physical significance. In the theory of phase transitions such a corner-stone was 
provided by the Ising model. It is natural to expect that a well-crafted dynamical version of the Ising model could 
play a key role in the development of non-equilibrium statistical mechanics. Driven lattice gases (Katz, Lebowitz and 
Spohn 1984; for a review see Schmittmann and Zia, 1995) provide such models; they represent particles hopping on a 
lattice and interacting through hard-core exclusion. The particles are subject to an external field that induces a bias 
in the hopping rates resulting in a stationary macroscopic current in the system. Due to this current, the microscopic 
detailed balance condition is violated and the system is generically in a non-equilibrium stationary state. 

The Asymmetric Simple Exclusion Process (ASEP) in one dimension is a driven lattice gas that can be viewed as 
a special case of the general class of models defined by Katz, Lebowitz and Spohn. This stochastic particle system 
was simultaneously introduced as a biophysical model for protein synthesis on RNA (MacDonald and Gibbs 1969) 
and as a purely mathematical tool for the study of interaction of Markov processes (Spitzer 1970, Liggett 1985, 
Spohn 1991, Liggett 1999). Subsequently, the ASEP has been used to study a wide range of physical phenomena : 
hopping conductivity in solid electrolytes (Richards 1977), transport of Macromolecules through thin vessels (Levitt 
1973), reptation of polymer in a gel (Widom et al. 1991), traffic flow (Schreckenberg and Wolf 1998), surface growth 
(Halpin-Healy and Zhang 1995, Krug 1997), sequence alignment (Bundschuh 2002) and molecular motors (Klumpp 
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and Lipowsky, 2003). 

From a theoretical point of view, the ASEP plays a fundamental role in the study of non-equilibrium processes 
(Krug 1991). Many exact results for the ASEP have been derived using two complementary approaches, the Matrix 
Product Ansatz and the Bethe Ansatz (for a review see Derrida 1998, Schiitz 2001), the relation between these two 
Ansatze being still a matter of investigation (Stinchcombe and Schiitz 1995; Alcaraz and Lazo 2004). 

The Matrix Product Ansatz, inspired from the quantum inverse scattering method (Faddeev 1984), is based on a 
representation of the components of the steady state wave function of the Markov operator in terms of a product of 
non-commuting operators. This Matrix Product technique, initially introduced for the ASEP with open boundaries 
(Derrida et al. 1993), has proved to be a very efficient tool to calculate steady state properties such as equal time 
correlations in the steady state, current fluctuations (Derrida et al. 1997), and large deviation functionals (Derrida 
et al. 2003) . This technique has been used to prove rigorously that the invariant distribution of the ASEP with two 
classes of particles is not a Gibbs measure (Speer 1993). 

The second approach, which consists in applying the Bethe Ansatz to a non-equilibrium process such as the ASEP, 
is due to D. Dhar (1987). The Markov matrix that encodes the stochastic dynamics of the ASEP can be rewritten in 
terms of Pauli matrices; in the absence of a driving field, the symmetric exclusion process can be mapped exactly into 
the Heisenberg spin chain. The asymmetry due to a non-zero external driving field breaks the left/right symmetry 
and the ASEP becomes equivalent to a non-Hermitian spin chain of the XXZ type with boundary terms that preserve 
the integrable character of the model. The ASEP can also be mapped into a six vertex model (Baxter 1982, Kandel 
et al., 1990). These mappings suggest the use of the Bethe Ansatz to derive spectral information about the evolution 
operator, such as the spectral gap (Gwa and Spohn 1992, Kim 1995, Golinelli and Mallick 2004a, 2005a) and large 
deviation functions (Derrida and Lebowitz 1998; Derrida and Appert 1999; Derrida and Evans 1999). 

The aim of the present work is to provide the reader with an introduction to integrability methods applied to the 
ASEP and to describe some of the exact results derived with the help of the Bethe Ansatz. Our presentation assumes 
little prior knowledge of these methods and we put emphasis on the algebraic properties of the model. In fact, the 
ASEP on a periodic ring is one of the most elementary systems to which the Bethe Ansatz can be applied : it is a 
discrete and classical model of interacting particles with the simplest possible interaction, hard-core exclusion, and 
with a conservation law (the total number of particles is constant). 

The layout of this work is as follows : in section |TT1 we describe the natural symmetries (by translation, reflection 
and charge-conjugaison) of the ASEP on a periodic ring; in section IIIII we review the solution of the ASEP by 
algebraic Bethe Ansatz and derive the Bethe equations that determine the spectrum of the model. In the case of the 
totally asymmetric simple exclusion process, the analysis of the Bethe equations can be carried out very precisely, even 
for finite size systems : in section IIVI we explain the procedure for solving the TASEP Bethe equations, determine 
the spectral gap, and show that a hidden symmetry of these equations allows to predict the existence of unexpected 
spectral degeneracies in the spectrum. In section |Vl we discuss the method for calculating the large deviation function 
of the total current; in particular, we show that the Gallavotti-Cohen relation manifests itself as a symmetry of the 
Bethe equations. In the last section IVTl we review the applications of integrability techniques to some variants of the 
ASEP. 

II. BASIC PROPERTIES OF THE EXCLUSION PROCESS 
A. Definition of the model : Tiie Markov matrix 

We consider the exclusion process on a periodic one dimensional lattice with L sites (sites i and L + i are identical). 
A lattice site cannot be occupied by more than one particle. The state of a site i {1 < i < L) is characterized by a 
Boolean number = or 1 according as the site i is empty or occupied. A configuration C can be represented by 
the sequence (ti, r2, . . . , tl). The system evolves with time according to the following stochastic rule (see Figure [1]): 
a particle on a site i at time t jumps, in the interval between t and t + dt, with probability p dt to the neighbouring 
site j -|- 1 if this site is empty {exclusion rule) and with probability q dt to the site z — 1 if this site is empty. The jump 
rates p and q are normalized such that p + q = 1. In the totally asymmetric exclusion process (TASEP) the jumps 
are totally biased in one direction (p = 1 or g = 1). 

We call Pt{C) the probability of configuration C at time t. As the exclusion process is a continuous-time Markov 
process, the time evolution of Pt{C) is determined by the master equation 

^R{C) = Y,M{C,C')PtiC') = J2 {Mo{C,C') + AhiCC) + M_,{C,C'))p,iC') . (1) 
c c 

The Markov matrix Al encodes the dynamics of the exclusion process: the non-diagonal element Mi(C,C') represents 
the transition rate from configuration C to C where a particle hops in the forward (i.e., anti-clockwise) direction, the 
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FIG. 1: The Asymmetric exclusion process on a periodic ring. 



non-diagonal element M_i(C,C') represents the transition rate from configuration C to C where a particle hops in 
the backward {i.e., clockwise) direction. The diagonal term Mq{C,C) = — X^c^^c + i^)) represents 

the exit rate from the configuration C. 

The matrix M is a real non-symmetric matrix and, therefore, its eigenvalues (and eigenvectors) are either real 
numbers or complex conjugate pairs. A right eigenvector P is associated with the eigenvalue E of M if 

MP = EP . (2) 

On a ring the total number of particles is conserved. For a given value of A^, the dynamics is ergodic (i.e., 
M is irreducible and aperiodic). The Perron- Frobenius theorem implies that is a non-degenerate eigenvalue and 
that all other eigenvalues have a strictly negative real part; the relaxation time of the corresponding eigenmode is 
T = — 1/Rc(i?). The right eigenvector corresponding to the eigenvalue i? = is the stationary state : for the ASEP on a 
ring the steady state is uniform and the stationary probability of any configuration C is given by P{C) = Nl{L~N)l/Ll. 

Remark : A configuration can also be characterized by the positions of the A^ particles on the ring, (xi, X2, . . . , xn) 
with 1 < xi < a;2 < ■ ■ • < xn < L. With this representation, the eigenvalue equation Q becomes 

EP{xi, . . .,xn) = 

^p[P{xi,...,Xi^l, Xi-1, Xi+i,...,XN)- Pixi,...,Xn)] + 

i 

'^q[P{xi, . . . ,Xj^i, Xj + 1, Xj+i,...,XN)- Pixi,...,XN)] , (3) 
j 

where the sum runs over the indexes i such that Xi-i < Xi — I and over the indexes j such that Xj + 1 < Xj+i ; these 
inequalities ensure that the corresponding jumps are allowed. 

The knowledge of the spectrum of the matrix M and of the associated eigenvectors provides a full description 
of the dynamics of the model. Such a spectral analysis is very similar to that carried out in Quantum Mechanics. 
One should however keep in mind that in the present case the evolution operator is non-hermitian : its eigenvalues 
are complex numbers and the associated eigenvectors do not form an orthogonal basis in the configuration space. 
Moreover, because M is not symmetric, its left eigenvectors are different from the right eigenvectors. 

B. Natural symmetries of ASEP on a periodic ring 

Before diagonalizing the Markov matrix M by Bethe Ansatz, we present some invariance properties of the ASEP, 
that can be described using elementary methods. The exclusion process on a ring displays symmetry properties under 
translation, right-left reflection, and particle/hole exchange (the last two symmetries play a role analogous to parity 
and charge conjugaison in quantum mechanics). These symmetries reveal some intrinsic formal properties of the 
Markov matrix, and are very useful to reduce the computation time when numerical diagonalization is performed. 
Finally, these natural symmetries will allow us to define some conserved quantities (analogous to quantum numbers) 
and to predict degeneracies in the spectrum. 
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1. Translation invariance 



The translation operator T shifts simultaneously all the particles one site forward: 

T{ti,T2,...,Tl) = {tl,Ti,T2,...,Tl-i). (4) 

Because of the periodicity we have = 1. Thus, the eigenvalues (or impulsions) fc of T are simply the L-th roots of 
unity: 

k = e^''^™/-^ for TO = 0, . . . , L - 1 . (5) 

We denote by T(fc) the eigenspace of T corresponding to the impulsion k. The projection operator 7^ over this 
eigenspace is given by 

L-l 

Tfc = - ^ k-^T^ . (6) 

The complex conjugation transforms a T-eigenvector of eigenvalue k into an eigenvector with eigenvalue k* = k~^, 
i.e., T{k*) = [T{k)]*. 

The ASEP on a periodic lattice is translation invariant, i.e., 

MT = TM. (7) 

The matrix M and the translation operator T can therefore be simultaneously diagonalized: M leaves each eigenspace 
T(fc) invariant. We denote by Sp{k) the set of the eigenvalues of M restricted to T(fc). Using complex conjugation 
we obtain the property 

Sp(r) = [Sp(A;)]* . (8) 



S. Right-Left Reflection 
The reflection operator R interchanges the right and the left and is defined by 

R{ti,T2,...,Tl) = {tl,TL-1,...,T2,Ti). (9) 

We have = 1; the eigenvalues of R are thus r = ±1. The reflection R reverses the translations, i.e., 

RT^T-^R, (10) 

and transforms a T-eigenvector of eigenvalue k into an eigenvector of eigenvalue k* . The operator R does not commute 
with the Markov matrix because the asymmetric jump rates (p ^ q) are not invariant under the exchange of right 
and left. Writing explicitly the dependence of M on the jump rates, we have 

RM^p,g)R-' = M(,,p) . (11) 

Thus, in general the reflection is not a symmetry of the ASEP (only the symmetric exclusion process is invariant 
under reflection). 



3. Charge conjugation 

The charge conjugation operator C exchanges particles and holes in the system, i.e., a conflguration with N particles 
is mapped into a configuration with L — N particles: 

C{ti,T2, . . . ,Ti) = (1 - Tl, 1 - T2, . . . , 1 - Ti) . (12) 

The operator C satisfies the relations: 

C2 = 1, CT = TC, CR = RC. (13) 

By charge conjugation, particles jumping forward are mapped into holes jumping forward. But holes jumping forward 
are equivalent to particles jumping backward. Writing explicitly the dependence of M on the number N of particles 
and on the jump rates, we thus have 

CM^N,p,g)C-^ = M^L-N,g,p) ■ (14) 

Wc notice that the number of particles is conserved by C only at half filling {L = 2N). Thus, the charge conjugation 
is not a symmetry of the ASEP except for the symmetric {p = q) exclusion process at half filling. 
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4- CR symmetry 

Using equations (fTTj) and p4|) to combine the charge conjugation C with the reflection i?, we obtain 

(Ci?)M(^,,p^,)(Ci?)-i = M^L-N,p.q) ■ (15) 

Thus, for a given L, the CR operator maps the ASEP with N particles into the ASEP with the same jumping rates 
but with L — N particles. This implies that the spectrum of M for N particles is identical with the spectrum for 
L — N particles (because CR transforms eigenvectors of M{^N^p^q) into eigenvectors of M(2^_jv,p,g))- 
For the ASEP model at half filling, i.e., L = 2N, the CR operator constitutes an exact symmetry : 

(CR) M = M (CR) . (16) 

The ASEP at half filling is therefore invariant under each of the two symmetries, translation T and CR. Note that 
they do not commute with each other; we rather obtain from equations (I10|l and ()13|) 

{CR)T = T-\CR). (17) 

Hence, the CR transformation maps the subspace T(fc) into T(fc*) and thus Sp(fc) — Sp(fc*) . For k ^ ±1, the 
subspaces T(fc) and T(fc*) are distinct and the corresponding eigenvalues of M are doubly degenerate. For k — ±1, 
the CR transformation leaves T(fc) invariant. Thus, the natural symmetries do not predict any degeneracies in the 
sets Sp(±l). However, each of the subspaces T(±l) is split into two smaller subspaces which are invariant under CR 
and on which CR = ±1. For the ASEP at half filling, we find from Equation ^ that the set Sp(fc) is self-conjugate, 
for all fc, and degenerate with Sp(fc*) 

Sp(fc) = [Sp(fc)]* = Sp(fc*) . (18) 
This means that Sp(fc), for all k, is made only of real numbers or of complex conjugate pairs. 



III. ALGEBRAIC BETHE ANSATZ FOR THE EXCLUSION PROCESS 



A. The algebra of local update operators 

The space of all possible configurations C is a 2^ dimensional vector space that we shall denote by S. Each site i 
is either occupied or empty : we shall represent its state by the local basis (|l)i, |0)i) of the two dimensional space 
Cf. The total configuration space S is thus given by 



S 



L y 



(19) 



The natural basis of the configuration space S is the tensor product (|l)i, |0))i (g) . . . (Xi (|1)l, |0))l- On a periodic 
lattice, the number N of particles is conserved by the dynamics. The total number of configurations for N particles 
on a ring with L sites is given hy Af — L\/[N\{L — A^)!]. The allowed configurations of the system span a subspace of 
S of dimension M . The Markov matrix can be expressed as a sum of local operators that update the bond (i, i + 1) : 



(20) 



i=l 



The update operator M^^i+i is given by 



/O 
-p g 
Q p -q Q 

Vo 



11+2 



(21) 



where Ij is the 2x2 identity matrix acting on the site number j. We emphasize that M^.i+i is a 2^ x 2^ matrix that 
acts trivially on all sites different from i and i + l and updates the bond (i, i + 1) according to the local dynamical 
rules of the exclusion process. We also define the local permutation operator Pi^i+i 



P, 



1,4+1 



10 
10 
10 
1 



Li+2 • ■ 



(22) 
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The local update operators satisfy the following relations: 







(23) 




= -M,,,+i 


(24) 




= pqM^^i+i 


(25) 




= pqMi^i+i 


(26) 


[Mi,i+i,Mjj+i] 


= for \i-j\>2. 


(27) 



These identities define a Temperley-Lieb algebra : this property of reaction-diffusion processes was emphasized by 
Alcaraz et al. (1994) and plays a key-role in the integrability of the ASEP. 



B. The ASEP eis a non-Hermitian spin chain 

The dynamics of the ASEP is entirely encoded in its Markov matrix M that governs the time evolution of the 
probability measure on the configuration space. We now show that the local operator Mj^j+i that updates the bond 
{i, i + located in the bulk of the system can be expressed with the help of the Pauli matrices. We recall that Pauli 
matrices are given by 

We also need the following operators 

After identifying the spin-i basis (| t), | i)) of with the local basis (|l)i, |0)j) of the two-dimensional space C? 
associated with the site i, we define an action of the Pauli matrices on the ASEP configuration space <S as follows : 

5," = li . . . lj_i O 6"" O li+i O . . . Iz, with a = x,y,z,+,- . (30) 

We observe that the local update operator Mj j+i can be written as 

M,,,+i = pS-St+, + qS+Sr,, + ^8^3^+, - ^^{3^ - 3^^,) - ^ , (31) 

where 1 represents the identity operator. The Markov matrix of the ASEP on the periodic lattice of size L is thus 
given by 



1=1 ^ 



i (32) 



(we recall that the site i -|- 1 is the same as the site number 1). The Markov matrix is therefore expressed as a spin 
chain on the periodic lattice. This spin chain is non-Hermitian when p q. For p = q = 1/2, the Markov matrix M 
is identical with the antifcrromagnctic Hciscnbcrg XXX spin chain which was exactly solved by Hans Betlic (1931). 
More generally, a similarity transformation allows us to map exactly the Markov matrix M on an antiferromagnetic 
XXY quantum spin chain with hermiticity breaking boundary terms. Following Essler and Rittenberg (1996), we 
consider the operators 

Z7i = li . . . ( Q J_i ] Ij+i . . . li, , (33) 



where we have defined 



The operator H defined as 



a = J-. (34) 



H={a + a-')il[uA(M + -l](Y[uA , (35) 
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is then given by 

^ = I E (s^Sf+i + SfSf+i + ^^-^f^S^SUi) + c^^SlSt + a-^S+S^ + ^±^SISI . (36) 

In this representation, we can formaUy rewrite H as the XXY spin chain on a periodic lattice with the following 
twisted boundary conditions 

Sl+i = SI. (37) 



C. The Yang-Baxter equation 

We shall now derive the Yang-Baxter equation which will allow us to define in the subsection IIII Di a one-parameter 
family of commuting operators that contains the Markov matrix M . Such a family of commuting operators ensures 
the existence of a sufficient number of conserved quantities that fully label the states of the ASEP. This property is 
the key to the integrability of the ASEP. 

Let us consider two auxiliary sites labeled as a and 6, each of them having two possible states (|l)a,|0)a) and 
(|l)b, |0)f,) respectively. The four possible states of a and h span a four dimensional complex vector space denoted as 
A® B. For any given value of the complex number A (called the spectral parameter) , we define an operator £afc(A) 
which acts on this tensor space as follows 

Cab{X) = -Pafc (1 + \Mab) , (38) 



where Mat and Pab represent the jump and the permutation operators between a and b as defined in equations (j2ip 
and respectively. We emphasize that the auxiliary sites a and b should not be viewed as neighbouring sites on a 
given lattice but rather as an 'abstract' pair of sites related by non-local jump and exchange operators Mat and Pat- 
We now consider three auxiliary sites a, b and c, and prove the Yang-Baxter relation : 

Cab{v)Ccb{>^)Cca{p) = Ccaip^CcbWCabiv) for V=- ^ , ^ . ■ (39) 

1 - ^ + pqX^J. 

We shall need the fact that the jump and permutation operators between the auxiliary sites satisfy the following 
relations analogous to 



we have 






p 

^ c,a 


Ma,b 


= M^^bPc^a 




Mlb 


= -Ma,b 


M^,a Ma.b 


Mc.a 


= pqMc.a 


Ma,b M^,a 


Ma.b 


= pqMa.b ■ 



(40) 
(41) 
(42) 
(43) 

In order to derive the Yang-Baxter equation, we first remark that 

CabHCcbWCM = Pab{l + yMab)Pcb{l + \M,b)Pca{l+ tlM,a) 
= Pab (1 + l^Mab) PcbPca (1 + AA^) (1 + M^ca) 

= PabPcbPca (1 + vM,a) (1 + AAfafc) (1 + t^M.a) ■ (44) 

Similarly, we have 

CMC,b[\)Cab[l^) = PcaPcbPab (1 + M^afc) (1 + XM^a) (1 + I^Mab) ■ (45) 

Now, with the help of the Temper ley-Lieb algebra, equations (|23ff27)) . we obtain 

(1 + l^M^a) (1 + XMab) (1 + /iMca) 
= 1 + (/i + u)Mca + XMab + IJivMl^ + XfiMabM^a + XvM^aMab + Xfil^M^aMabMca 

= + I^LV+ pqXfliy)Mca + XMab + XflMabMca + XvMcaMab ■ (46) 

Similarly we have 

(1 + ^lMab) (1 + XM,a) (1 + vMab) 
= l + {^i + V-pLV+ pqXtJ.l^)Mab + XM^a + XfiMabM^a + XvM^aMab ■ (47) 

Using the fact that PabPcbPca = PcaPcbPab wc complete the proof of equation 
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D. The Transfer Matrix 



We are now ready to apply the Algebraic Bethe Ansatz to the ASEP (for introduction to this subject, see Faddeev 
1984, de Vega 1989, and Nepomechie 1999). We introduce an auxiliary site (denoted as site 0) which can be in two 
states |l)o or |0)o- These two states span a two dimensional vector space. A, the auxihary space. We define, as in 
equation pS)) . a local transfer operator Cio{X) between the site i of the ASEP ring and the auxiliary site 0. This 
operator, that we shall denote as Ci(\), can be represented as a 2 x 2 operator on the vector space A, 

A(A) = ( 'ti'l ) , (48) 



c(A) d{\) 

where the matrix elements a(A), 6(A), c(A) and d{X) are themselves 2^ x 2^ operators that act on the configuration 
space S. These operators act trivially on all sites different from i and are given by 

a(A) = li (8) . . . (g) J gA ) ® -^'+1 ■ ■ • ' ^^^^ 
b{X) = U(g)... (g) ^ [J ^ ® . . . 1l , (50) 

c(A) = li ^ . . . l,_i (g Q ^ 0*^^ ) ^^+1 • ■ ■ ' 

d(A) = li ® . . . <E,(P^ 5 ] ® (52) 



We now consider the operator 



n{iy) = l + l^Ma>,a, (53) 



that acts on AiSi A' where the auxiliary spaces A and A' are the configuration spaces of the auxiliary sites a and a' . 
In the basis (|la, la'), |la, 0^'), |0a, 1^'), |0a, Oa')) of A(^A', the operator TZ{i^) is represented by a 4 x 4 scalar matrix : 

10 

1 

Using equation (|39|). we find that the operators £i(A) satisfy the Yang-Baxter equation : 

7^(^y) [£, (A) (g, C\ (p)] = [£, (m) (A)] Tliv) with = ^ , (55) 

1 — ^ + pqA^ 

where /!i and C i are interpreted as 2 x 2 matrices acting, respectively, on A and A' with elements that are themselves 
operators on S. Their tensor product is thus a 4 x 4 matrix, acting on A® A! with matrix elements that are operators 
on S. 

The monodromy matrix is defined as 

T{\)=Ci{\)C2{\)...Cl{\). (56) 

where the product of the CiS has to be understood as a product of 2 x 2 matrices acting on A with non-commutative 
elements. The monodromy matrix T(A) can thus be written as 



_ / A{\) B{X) 
^^^>-\ C(A) D{X) 



(57) 



where A, B, C and D are operators on the configuration space S. By applying equation (f55|) to all the lattice sites i, 
we find that the monodromy matrix satisfies a similar relation 

7e(z/) [r(A) (g T(^)] = [Tifi) (g ^(A)]7^(^/) with = ^ , (58) 

1 — fi + pqXfi 
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Taking the trace of the monodromy matrix over the auxiUary space A, we obtain a one-parameter family of transfer 
matrices acting on S 

t{\)=Ti-A{T{X))=A{X)+D{X). (59) 

By taking the trace of equation ((58|) over A(E) A' and using the fact that TZ{i') is generically an invertible matrix, we 
deduce that the operators t{X) form a family of commuting operators (see e.g., Nepomechie, 1999). In particular, this 
family contains the translation operator T — t{0) (that shifts all the particles simultaneously one site forward) and 
the Markov matrix M = t'{0)/t{0). 

The equation when written explicitly, leads to 16 quadratic relations between the operators ^(A), B{X), C(A), 
D{X) and the operators A(/i), B{ii), C{ii), D{ii). In particular, we have 

c{x)c{ii) = c{^l)c{x), (60) 
A{x)C{^i) = {\-qv)A{^Ji)C{x)+pvC{^i)A{x), (6i) 

D{X)C{^l) = piyC{fi)D{X) + (1 - pu)D{^i)C{X) . (62) 
These relations are used in the next subsection to construct the eigenvectors of the family of transfer matrices t{X). 



E. The Bethe equations 



Using Algebraic Bethe Ansatz, the common eigenvectors of this family are explicitly constructed by actions of the 
C operators on the reference state ft, defined as 

= |0i02...0l). (63) 

The state ft corresponds to a configuration where all the sites are empty. The operators A{X) and D{X) have a simple 
action on ft, namely 

A(-]n^X^n and D(-]n^n. (64) 



Now, for any N < L, we define the vector 

\zuZ2...z^)=C^^y..c[^y[^^ n, (65) 

where Zi, Z2 ■ ■ ■ are complex numbers. Each operator C creates a particle in the system, and the state \zi, Z2 ■ ■ ■ Zn) 
is a linear combination of configurations with exactly N particles. This vector is an eigenvector of the operator t{X/q) 
(for all values of A) and in particular of the Markov matrix M, provided the pseudo-momenta zi, Z2, . . ., zjy satisfy 
the Bethe equations (Gwa and Spohn 1992) : 

N _ 

= (-1)^-1 TT ^^'^' ' for 1 = 1. ..N. (66) 
f-J^ pziZi - Zi+q 

These equations are derived by considering the action of the operators t(X/q) on |zi, Z2 ■ ■ ■ zn) : 

Kt)'— >^^e)Kf)-K7)<7)-K^)Kf)-K7)<7)- 

Applying repeatedly equations (|6ip and (|62p to the first term and to the second term on the r.h.s. of this equation, 
respectively, we obtain a linear combination of \zi, Z2 ■ ■ ■ zn) and of N vectors |zi, . . . z^-i. A, Zi+iz^r) for i = 1, . . . , N . 
If we want |zi, Z2 ■ ■ ■ zm) to be an eigenvector of t{X/q), the coefficients of these N vectors must vanish . The vanishing 
conditions for the coefficients of the unwanted terms provide a set of N relations that must be satisfied by the complex 
numbers zi, Z2 ■ ■ ■ zn- These relations, when written explicitly, are precisely the Bethe equations (|66p . 
The corresponding eigenvalue of t{X/q) is given by 

Pf\ \ rj P>^z,-X + q pXz,- z, + q 

£{X]zi,z2. ■ -zn) = YY — hA 11 — . (68) 
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Using the Bethe equations ([66)) . we find that £{X) is a polynomial in A of degree L — N. The eigenvalue of the Markov 
matrix is given by 



N N 



i=l i=l 

The Bethe equations (|66l) for the ASEP were first derived and analyzed by Gwa and Spohn (1992), who used the 
coordinate Bethe Ansatz, a method more elementary and quicker than the one presented here. We have chosen to 
describe here the more sophisticated Algebraic Bethe Ansatz technique, because, in our opinion, the ASEP provides 
the simplest pedagogical example of this method. Besides, the Algebraic Bethe Ansatz becomes unavoidable to study 
exclusion processes with different classes of particles and systems with open boundaries (de Gier and Essler 2005). 
Finally, the Algebraic Bethe Ansatz allows to construct a Matrix Product Ansatz involving quadratic algebras and to 
draw a rigorous relation between the two complementary approaches used to analyze the ASEP (Golinelli and Mallick 
2006). 

IV. ANALYSIS OF THE BETHE EQUATIONS OF THE TASEP 

In the limit of a totally asymmetric exclusion process the Bethe equations become particularly simple and can be 
analyzed thoroughly. In this section we present some exact results that are valid for the TASEP. Substituting p = 
and q = I in the Bethe equations (j66p . we obtain 

N 

= (-1^-^11131^ ^ = 1...A^. (70) 

j=l ^3 

In terms of the variables Zi — 2/zi — 1, these equations become (Gwa and Spohn 1992) 

(l-Z,)^(l + Z,)^-^ = -2^n|^ ^ = l,...,iV. (71) 

j=i 3 

We note that the right-hand side of these equations is independent of the index i: this property is true only for 
the TASEP where the Bethe equations decouple and can be reduced to an effective one-variable equation. The 
corresponding eigenvalue E of the Markov matrix M is given by 

N 

2E ^ ~N + Zj . (72) 

For the TASEP, it is possible to show that the associated Bethe eigenvector Pe is a determinant (Golinelli and Mallick 
2005b) 

PE{xu...,XN)^det{R), (73) 

where R is a N x N matrix with elements 

Rii,j) = j-^^ forl<z,j<iV, (74) 
(i — Ziy 

where (zi, . . . , zn) are the roots of the Bethe equations ([70|) . Similar formulae with determinants also appear in the 
expression of Pt calculated for an arbitrary initial condition Pq on an infinite open lattice (Schiitz 1997, Sasamoto 
and Wadati 1998, Rakos and Schiitz 2005) and on a periodic lattice (Priezzhev 2003). 

A. Layout of the solutions of the Bethe equations 

The solutions of the Bethe equations ([7T|) are the roots of the polynomial equation of degree L 

(1-Z)^(1 + Z)^-^ =y, (75) 
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where Y must be determined self-consistently by the r.h.s. of equation ((7T|) . In this subsection, we describe the 
geometrical layout of the roots of the polynomial (j75p for an arbitrary value of the complex number Y. 

The L solutions {Zi, . ..,Zl) of this equation display interesting geometric properties. We explain how to label 
these solutions so that each Z^iY) becomes an analytic function of the parameter Y in the complex plane with a 
branch cut along the real semi-axis [0, -|-cxd). 



-N+1 



-N 



N-1 



-L-1 



FIG. 2: Labelling the roots of the Bethe equations. Here L — 15, N = 6, tf) = n/2 and r/rc — 0.8,1,1.2 (see text). The 
continuous curves are Cassini ovals. When r is fixed and (j) varies from to 2tt, each Zu slips counterclockwise along a part of 
the Cassini oval. When 4> is fixed and r varies from oo to 0, each Zk travels along a dashed curve from oo to points +1 or -1. 



A non-zero complex number Y can be written in a unique way as 

Y = r^ e'^ with < < 27r , 



(76) 



r being a positive real number. We emphasize that the argument has a branch cut along [0, +oo) : (j) has a discontinuity 
of 27r when Y crosses the positive real axis. For a given value of r, the complex numbers Zk belong to the generalised 
Cassini oval defined by 



mz + i 



(77) 



where p = N/L is the filling of the system. As shown in Fig. [21 the topology of the Cassini ovals depends on the value 
of r with a critical value 



(78) 



• for r < rc, the curve consists of two disjoint ovals with N solutions on the oval surrounding +1 and L — N 
solutions on the oval surrounding —1. 

• for r = rc, the curve is a deformed Bcrnouilli Icmniscatc with a double point at Zc ~ 1 — 2p. 

• for r > rc, the curve is a single loop with L solutions. 

The Cassini ovals are symmetrical only if p = 1/2. 

In order to label the solutions, we start by considering the limit r — > oo for a given (f>. Equation (|75p then becomes 



exp[i((/) - Ntt)] 



The solutions Zk are given by 



r exp 



NTT + 2(k- 1)tt] 



with k = I, . . . , L . 



In other words, the Zk are regularly distributed along a large circle of radius r with 

2n. 



<- Ntt (f>- Nn 

— < arg Zl < ■ ■ ■ < arg Zl < 



L 



(79) 



(80) 



(81) 



12 



This labelling, obtained for large r, is extended by analytic continuation to all values of r, keeping 4> fixed. The loci 
of the Zfc are plotted in Fig. [5] (dashed curves): they are orthogonal to the Cassini ovals. A singularity appears along 
the branch (p — Q because Zi and Zm+i collapse into each other at the double point Zc when r = Tc] we circumvent 
it by choosing = 0+. 

With this labelling, the solutions are ordered along the Cassini ovals. Moreover, when r < Tc, the solutions 
{Zi, . . . , Zn) group together on the right oval and (Zjv+i, . . . , Zl) on the left oval. 



B. Procedure for solving the Bethe Equations 

The special 'one-body' structure of the TASEP Bethe equations leads to the following self-consistent method for 
solving them. In the following sections, we shall show how this procedure can effectively be used to calculate the 
spectral gap of the model and to determine some arithmetical global properties of the spectrum. 

• SOLVE, for any given value of Y, the polynomial equation (1 — Zi)^ {1 + Zi)^^^ = Y . The L roots 
Zi{Y) . . . Zl{Y) of this equation are located on a Cassini Oval 

CHOOSE N roots Zc(i)(F),... Zc(jv)(5^) amongst the L available roots, i.e., select a choice set c : 
{c(l),...,c(iV)}c{l,...,L}. 



• SOLVE the self- consistent equation Ac{Y) = Y where 

N 



A,(r) = -2^n 



ZcU){Y)-l 



The self-consistent equation Ac{Y) = y is a fixed point equation. It is solved in practice by choosing iteratively 
a new Y' given by Y' = Ac{Y) and going back to the first step until convergence is reached. 

• CALCULATE from the fixed point value of Y, the Zc(j)(F)'s and the energy corresponding to the choice set 



c : 



N 



2E,(Y) = -N + \^Z,(,){Y). 

Consider the choice function co(j) — j that selects the N fugacities Zi with the largest real parts. The ^-function 
and eigenvalue associated with this choice function are given by 

Ao{Y) = -2^l[^—, (82) 

N 

2Eo = -N + J2Zj. (83) 
i=i 

In (Golinelli and Mallick 2004a, 2005a), we derived the following exact combinatorial formulae for Ao{Y) and Eo{Y) 
for any finite values of L and N : 

fc=i ^ ^ 

^o(n = -f:(tN-l)^- (85) 

k=l ^ ^ 

These expressions are derived as follows: when F ^ 0, the TV roots of the equation ([75)1 with the largest real parts, 
i.e., [Z\, . . . , Z-^) converge to -1-1 (whereas the L — N other roots go to -1). Consider a contour 7, positively oriented, 
that encircles -1-1 such that for Y small enough (Zi, . . . , Ziq) are inside 7 whereas {Z^+i, ■ . . , Zl) are outside 7. Let 
h{Z) be a function, analytic in a domain that contains the contour 7; from the residue theorem, we deduce: 

^,,rr^ 1 r {\- Z)^-^{1 + Z)^-^-^{L-2N - LZ) 



^ ^(^^-^ - % (i-W-^z)^-^-y ■ 
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The functions log{Ao{Y)/Y) and Eo(Y) are expressed as such contour integrals by choosing h{Z) = log(ii^) and 
h{Z) — Z — 1, respectively. Equations ([M]) and ([55)1 result from expanding for small values of Y the denominator in 
equation thanks to the formula {A - Y)-^ = J2T=o Y'^/A'' valid for |r| < \A\. 

The expressions (j84p and (j85p are then analytically continued in C — [1, oo). In the limit L oo and with p fixed, 
we obtain from the Stirling formula 

MY) 1 r: 



where Tc was defined in equation (j78p and the polylogarithm function Li^ of index s is given by 

The function Lig is defined by the first equality on the whole complex plane with a branch cut along the real semi-axis 
[l,-|-oo); the second equality is valid only for \z\ < 1. Similarly, when L —^ oo, Eq(Y) can be expressed in terms of 
the polylogarithm function Li5/2- 

From equation ((84|) . we observe that the equation AgiY) = Y has the solution Y — that yields Zj = 1 for j < N. 
This solution leads to i?o = 0. The choice function co(j) = j thus provides the ground state of the Markov matrix. 
This fact has been verified numerically for small size systems (Gwa and Spohn 1992, Golinelli and Mallick 2004a). 

C. Calculation of the Gap 

The spectral gap, given by the first excited eigenvalue, corresponds to the choice ci{j) — j for j = 1, . . . , — 1 and 
ci{N) = + 1 (Gwa and Spohn 1992). The associated self-consistency function and eigenvalue are given by 

MY) = Ao(r)|^|^, (89) 
2Ei = 2Eo + {Zn+i- Zn). (90) 
In the large L limit, Bethe equations for the gap become at the leading order 

Li3/2(-e"'') = 2iTT \{u + ifl'^ -{u- i)^/^] , where u is defined as ^ = -e"'' . (91) 

The solution of this equation is given by: 

u= 1.119068802804474... (92) 
This leads to the eigenvalue corresponding to the first excited state : 

g.^-27^ "™--- ±g^ig^, (93) 

We observe that the first excited state consists in pair of conjugate complex numbers when p is different from 1/2. 
The real part of Ei describes the relaxation towards the stationary state. The corresponding relaxation time scales 
as the size L of the system raised to the power 3/2; the dynamical exponent z of the ASEP is thus given by z = 3/2. 
This value agrees with the dynamical exponent of the one-dimensional Kardar-Parisi-Zhang equation that belongs 
to the same universality class as ASEP (Halpin-Healy and Zhang 1995). The imaginary part of Ei represents the 
relaxation oscillations and scales as L~^; these oscillations correspond to a kinematic wave that propagates with the 
group velocity 2p — 1 (Majumdar, private communication). The method described here can be extended to calculate 
the higher excitations of the spectrum above the ground state by considering other choice sets c. However, the 
precise correspondence between the choice set and the level of the excitation has not been studied at present, to our 
knowledge. 
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L 


iV 


m(l) 


m(2) 


m(6) 


m(20) 


m(70) 


2 


1 


2 










4 


2 


4 


1 








6 


3 


8 


6 








8 


4 


16 


24 


1 






10 


5 


32 


80 


10 






12 


6 


64 


240 


60 


1 




14 


7 


128 


672 


280 


14 




16 


8 


256 


1792 


1120 


112 


1 


18 


9 


512 


4608 


4032 


672 


18 



TABLE I: Spectral degeneracies in the TASEP at filling p = 1/2; L is the size of the lattice, N the number of particles, m{d) 
the number of multiplets with degeneracy d. 



p 


L 


iV 


m(l) 


m(2) 


m(3) 


m(4) 


m(5) 


m(15) 


1/3 


9 


3 


81 




1 










12 


4 


459 




12 










15 


5 


2673 




90 


15 








18 


6 


15849 




540 


270 




1 




21 


7 


95175 




2835 


2835 


189 


21 


1/4 


16 


4 


1816 






1 








20 


5 


15424 






20 








24 


6 


133456 






240 


36 




1/5 


25 


5 


53125 








1 




2/5 


15 


6 


4975 


15 











TABLE IL Examples of spectral degeneracies in the TASEP at filling p ^ 1/2; m{d) is the number of multiplets with degeneracy 
d. 

D. Spectral degeneracies due to a hidden symmetry of the Bethe equations 

Although the steady state of the ASEP and the lowest excitations have been extensively studied, little effort has 
been devoted to investigate global spectral properties of the Markov matrix. In this section, we show that the 
Bethe equations of the TASEP possess an invariance property under exchange of roots that implies the existence 
of unexpected multiplets in the spectrum of the TASEP evolution operator. This hidden symmetry of the Bethe 
equations allows to predict combinatorial formulae for the orders of degeneracies and the number of multiplets of a 
given order of degeneracy. 

A numerical diagonalization of the Markov matrix reveals a striking global property of the TASEP spectrum : the 
existence of large degeneracies of very special orders. Degeneracies of orders 2, 6, 7, 20... appear in the TASEP 
spectrum at half- filling (Table H]). This property does not follow from the natural symmetries of the exclusion process 
(discussed in section Hi Bp which suggest that the spectrum should be composed only of singlets for impulsion = ±1 
and of singlets and doublets for k ^ ±1. Degeneracies of higher order have no reason to appear and therefore should 
not exist generically; however, multiplets do exist in spectrum of the TASEP at half-filling (Table H]) as well as in the 
TASEP at arbitrary filling (Table |TT|. 

These spectral degeneracies have an arithmetical origin and appear only when L (the number of sites) and N (the 
number of particles) are not relatively prime. We thus define 

S = gcd(L, N) . (94) 

The L Bethe roots of the polynomial l|75|) can be grouped into S disjoint packages, each of cardinality L/S. The roots 
composing the package Vs have the indices {s,s + S,s + 26,...,s + L — 6} with 1 < s < 6. Consider a choice set c {i.e., 
as explained in section HVBi a choice of N roots amongst the L available ones). Suppose there exist two packages Vs 
and Vt such that 

VsCc and n c = . 

Consider the choice set c = {c\Vs)^Vt obtained by replacing in the choice set c the roots composing Vs with the roots 
Vt- It can be proved (Golinelli and Mallick 2005b) that this choice set c (obtained from c by exchanging Vs and Vt) 
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corresponds to the same self-consistent equation and to the same eigenvalue E as c. This property therefore defines 

Equivalence Classes amongst the ^ ^ ^ possible choice sets : two choice sets are equivalent if they are obtained from 

each other by 'Package-swapping'. 

Figure [3] illustrates the case of a system with L — \{) sites and N = 5 particles, where 5 Bethe roots must be chosen 
amongst the 10 solutions {Zi, . . . Zio) lying on a Cassini oval of a polynomial of degree 10. Here 6 = gcd(L, N) — 5, 
and therefore the roots can be grouped into five different packages : 

Vi - {Zu Ze} , V2 - {Z2, Zr} , V3 = {Z3, Zg} , V4 - {Z4, Zq} , V5 = {Z5, Z^} ■ 

These five packages correspond to the five rectangular boxes in figure [Sj each package contains two roots and the 
selected root is shown as black disk. In this figure, we give the example of a choice function c that selects the roots 
Zi, Zii Z-i^ Zg, Z4. Then, the choice function c that selects the roots Zi, Zg, Z4, Z5, Zio, is such that c = {c\V2) U V5, 
i.e., the two choice sets c and c are obtained from each other by exchanging the packages V2 and V5. They are thus 
equivalent in the sense defined above and correspond to the same eigenvalue E = l/2(— 5 + Zi + Z^ + Z4). 



26 



PI 



P2 



P3 



P4 



P5 



• 


Z2 


• 


Z3 





Z4 


• 


Z5 










• 




• 


Z9 
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C 



O 



o 
o 



o 



o 



C AND C CORRESPOND TO THE SAME EIGENVALUE 
FIG. 3: Example of two equivalent choice functions of the Bethe roots for the TASEP with L = 10 and N = 5. 

The calculation of the degeneracies tables is therefore reduced to a pure problem in combinatorics : given a choice 
function c, one has to count the number of choice functions that are equivalent to c, and the total number of equivalence 
classes. Indeed, if we suppose that there is a 'one to one' bijection between choice sets and solutions of the Bethe 
Equations (recall that the number J\f of possible choice sets is the same as the dimension of the matrix M), then the 
following correspondences hold : 

• equivalence classes by 'package swapping' ^ multiplets in spectrum 

• cardinality of a class ^ order of the multiplet 

• number of classes of cardinality d ^ number of multiplets of degeneracy d. 

In the half- filling case, the numbers in Table U are given by simple formulae (Golinelli and Mallick 2004b). The 
degeneracy order d and the number m{d) of multiplets of a given degeneracy d can be expressed as a function of a 
single integer r 



2r 



, m{dr) 




(95) 



where r takes all integral values in the range < r < . arbitrary filling, expressions for m and d have also been 
found but they depend on {L/6— 1) different parameters (Golinelli and Mallick 2005b). 
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V. CALCULATION OF LARGE DEVIATION FUNCTIONS 



The Bethe Ansatz has also led to an exact calculation of all the cumulants of the total current in the TASEP on 
a ring of size L with TV particles (Derrida and Lebowitz 1998, Derrida and Appert 1999). As explained below, the 
generating function of these cumulants can be expressed as the ground state of a suitable deformation of the Markov 
matrix. From this generating function, the large deviation function of the time averaged current is obtained exactly. 



A. A generalized master equation 

We call Yt the total distance covered by all the particles between time and time t and define Pt{C,Y) the joint 
probability of being at time t in the configuration C and having Yt — Y. A master equation, analogous to equation ([1]), 
can be written for Pt{C, Y) as follows : 

^P*(c,r) = |J](Mo(c,c')Pt(c',r) + Mi(c,c')Pt(c',y-i) + A/_i(c,c')Pt(c',r + i)) . (96) 

In terms of the generating function Ft{C) defined as 

oo 

i^t(C) = ^e-^^PtlCr), (97) 

Y=0 

the master equation (j96p takes the simpler form : 



^FtiC) - J2 (Mo(C,C') +eWi(C,C') +e-'^Af_i(C,C'))i^t(C') = ^ Af(7)(C,C')Ft(C') . 
c c 

This equation is similar to the original Markov equation ([1]) for the probability distribution Pt{C) but where the 
original Markov matrix M is deformed into M^j) which is given by 

M(7) = Mq + e^Mi + e-''M_i . (99) 

We emphasize that M(7), that governs the evolution of Ft{C), is not a Markov matrix for 7 7^ (the sum of the 
elements in a given column does not vanish). In the long time limit, t — > 00, the behaviour of Ft(C) is dominated by 
the largest eigenvalue A(7) of the matrix M{j) : 

Fi(C)^e^(^)*(C|A(7)), (100) 
where the ket |A(7)) is the eigenvector corresponding to the largest eigenvalue. Therefore, when t — > 00, we obtain 

--'^'^ ^J^^tiC) ^ e^'^''^^ (101) 



e' 



More precisely, we have 



lim ilog(eT'*'*) = A(7). (102) 

t^oo t 



The function A(7) contains the complete information about the cumulants of the total current Yt in the long time 
limit. The large deviation function G{j) for the total current j is defined as 



1 



Yt 



G(j) = lim - log Prob {^=j] ■ (103) 
The relation between A(7) and G{j) is derived as follows 

( e-^^' ) ~ e^(^)* ~ / e*(«(^)+^^Mj . (104) 
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By saddle-point approximation, we deduce that A(7) is the Legendre transform of G{j) 

A(7) = max (gU) + 7j) • (105) 

This relation allows one to obtain the large deviation function G(j) in the following parametric form, once A(7) is 
known : 

J = ^ and G(j) = A(7) - 7^ • (106) 
d7 d7 

The largest eigenvalue A(7) of the deformed matrix M{'j) is calculated by Bethe Ansatz. The Bethe equations now 
read 



N 

'■I 



(-l)"-^n '?'''^''t''^^ for ^ = l...iV, (107) 



and the corresponding eigenvalue of M{^) is given by 

N N 



N I 

E{r, zi,Z2...zn)^ pe^ 9^"'' Y.--^- (108) 



i=l 



B. The Gallavotti-Cohen symmetry 

We remark that the equations (|107|) and (|108p are invariant under the transformation 

1 

z 

7 — > 70 — 7 with 7o = log — . (109) 

P 

This symmetry implies that the spectrum of -^(7) and that of M(7o — 7) are identical. This functional identity is in 
particular satisfied by the largest eigenvalue of M and we have 

A(7) = A(7o - 7) • (110) 
Using equation (jlOSp . we deduce that this identity implies the following symmetry for the large deviation function 

GU) = G{-j) - 70J . (Ill) 

This relation is a special case of the general Fluctuation Theorem valid for systems far from equilibrium (Evans, Cohen 
and Morriss 1993; Evans and Searles 1994, Gallavotti and Cohen 1995). This relation, which manifests itself here 
as a symmetry of the Bethe equations, does not depend on the integrability of the system and can be derived from 
equation directly (Kurchan 1998, Lebowitz and Spohn 1999). 



C. The large-deviation function for the TASEP 

In the TASEP case, the Bethe equations (|107p are simpler and can be solved by using the procedure outlined 
in section IIVBI An exact power series expansion in 7 for A(7) is obtained by eliminating the parameter Y in the 
following two equalities : 



fe=l 



fe=l 



kN kL-\ 
kN I kL 
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Equation (|113p is inverted by writing Y formally as a power series in 7 and solving recursively for the coefficients up 
to any desired order; the function A(7) is then obtained by substituting Y{'-f) in equation (jll2p . 
These expressions allow to calculate the cumulants of Yt , for example : 



lim 

t — ^00 



limi^ 

t^oo t 



dA(7) 



d7 
d2A(7) 



d72 



_ N{L - N) 
7=0 L ~ 1 ' 

N'^ {2L - i)\ [N - [L - N)9 
7=0 ^ {L - 1)!2 {2N - 1)! (2L -2N - 1)! 



In the limit of a large system size, L — > 00, and with a fixed density p = L/N, these expressions reduce to 



lim 



(Yt) 



Lp{l~p), and lim _ 

t^oc t I 



M^.^L3/2[p(l-p)]3/2 



t^oo t 

Using the Stirling formula, equations (|112[) and (|113[) take the scaling form (Dcrrida and Lebowitz 1998) 



A(7) - 7ip(l -P)^\\ ^^IT^'f (7 V2^P(1 - 



(114) 
(115) 

(116) 

(117) 



where the expression for the scaling function $ does not depend on the parameters (L and p) of the model. From 
equation (jll7p . it can be shown that when |j — Lp[\ — p)\ <C i, the large deviation function G(j) can be written as 



with 



G(j) ^ 



p{l-p)jj(]-Lp{l^p) 



7rL3 



.(]__±i 
V oil 



P(l - P) 



2V3 
4\/^ 



for 

y|3/2 for 



y 



(118) 

(119) 
(120) 



This distribution is skew: it decays as a power law with an exponent 5/2 for y +cx) and with an exponent 3/2 for 
y -00. 



VI. APPLICATIONS TO RELATED MODELS 



In the previous sections, we have put emphasis on the Bethe Ansatz solution of the totally asymmetric exclusion 
process on a ring. In this special case, the Bethe equations have a particularly simple form and their analysis can 
be carried out thoroughly. In this last section, we briefly present results that have been obtained with the Bethe 
Ansatz for the partially asymmetric exclusion process (PASEP), for the exclusion process with an impurity, and for 
the exclusion process with open boundaries. 



A. The partially asymmetric exclusion process 

The Bethe equations ()66p for general hopping rates p and q can not be reduced to an effective one- variable polynomial 
equation as in the TASEP case and their analysis becomes much more involved. In particular, combinatorial expansions 
for the energy such as ([55]) valid for finite values of L and N are not known. However, a very thorough study has 
been carried out, in the limit L — > 00, by Doochul Kim and his collaborators in a series of papers (Noh and Kim 1994, 
Kim 1995, Kim 1997, Lee and Kim 1999; see also related works on the non-Hermitian XXZ chain by Albertini et al. 
1996, 1997). They developed a perturbative scheme that enabled them to calculate the finite size corrections of the 
gap and the low lying excitations of the asymmetric XXZ chain. They derived the crossover scaling functions of the 
energy gaps from the symmetric case p = q = 1/2 to the asymmetric case. In the continuous limit, these functions 
describe a crossover from the Kardar-Parisi-Zhang universality class to the Edwards- Wilkinson class. 

The same technique leads to the calculation of the large deviation function of the PASEP in the scaling limit (Lee 
and Kim 1999). The large deviation function takes the form 
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where the function H is the same as that of equation (|118p . The only modification that occurs for the PASEP is 
the rescahng factor p ~ q. Non-trivial differences appear in the subleading terms, but these corrections are likely 
to be model-specific and non-universal. Using this expression of the large deviation function, it is possible to verify 
explicitly that the Gallavotti-Cohen relation (jllip is satisfied (Lebowitz and Spohn 1999). However, an exact power 
series expansion of A(7) in terms of the deformation parameter 7 is not known. Such an expansion would lead to exact 
formulae for the cumulants of the current valid for any values of L and N. We believe, nevertheless, that expressions 
analogous to equations (|112l and I113P should exist for the PASEP because exact combinatorial formulae for the mean 
value of the current and its variance have been calculated by using the Matrix Product method (Derrida and Mallick 
1997). 

B. ASEP with an impurity 

The presence of a defective particle (say, an impurity) in the ASEP on a ring can generate a shock dynamically 
in the stationary state (Mallick, 1996). A trafhc-flow picture nicely illustrates this fact : if particles represent cars 
and the impurity a truck (which moves at a slower speed and is difficult to overtake), then the shock corresponds 
to a traffic jam. The phase transition to a shock can also be interpreted as a real space version of Bose-Einstein 
condensation (Evans 1996). 

The model with a single defect particle is drawn in Figure 3) The system contains N normal particles (denoted by 
1) and one impurity (denoted by 2). Each site is either occupied by a particle or by the impurity, or it is empty. The 
stochastic dynamical rules that govern the evolution of the system during the infinitesimal time step dt are 

10 —> 01 with rate 1 

20 02 with rate a (122) 
12 21 with rate f3 . 

All other transitions are forbidden. The phase diagram of this model can be determined exactly by the Matrix 
product method and consists of four main phases (see Figure |4]). In the two massive phases, the perturbation due to 
the impurity is of finite range and all correlations decay exponentially. In the non-massive phase, the impurity has a 
long range effect : correlations decay algebraically. In the shock phase, the presence of the unique impurity induces 
a phase separation between a low density region and a high density region. 

By using the Bethe Ansatz, Derrida and Evans (1999) have calculated the large deviation function of the total 
displacement of the defect particle. Their result leads to analytical formulae for the velocity and the diffusion constant 
of the impurity in the various phases. In particular, in the shock phase, the location of the defect can be identified with 
the position of the shock and the large deviation function provides complete information on its statistical behaviour. 




FIG. 4: The Asymmetric exclusion process on a periodic ring with an impurity. 



A more direct way to produce a shock in a system is to introduce a slow bond that particles cross with a rate r < 1 
(Janowski and Lebowitz 1992, 1994). For r less than a critical value that depends on the mean density of particles, the 
system exhibits a phase separation. This model, in spite of many efforts, has not been solved exactly and is probably 
not an integrable system (it does not satisfy any known criterion for integrability) . However, G. Schiitz (1993) has 
solved, using Bethe Ansatz, a variant of this model in which particles move deterministically on a ring of length L 
(with L even) with a single defect across which the particles jump with probability r < 1. In this dynamics, the time 
evolution consists of two intermediate time steps in which even sites and odd sites are updated simultaneously. 
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C. The open boundary model 

The asymmetric exclusion process with open boundaries can be viewed as a model for a driven lattice gas in contact 
with two reservoirs. In the bulk, the dynamics of the ASEP is the same as that shown in Figure 1; at the boundaries 
particles can enter or leave the system with various input/output rates. The totally asymmetric version of this model 
is illustrated in Figure [5] : particles in the bulk jump to to right with rate 1; a particle can be added with rate a if 
the site 1 is empty; a particle can be removed with rate /3 if the site L is occupied. 

The open ASEP undergoes boundary induced phase transitions (Krug 1991). The density profile and the current 
can be calculated exactly by using a Matrix product representation for the steady state (Derrida et al. 1993; see 
also Schiitz and Domany 1993 for an alternative derivation). In the large system size limit, the expressions for the 
current and the density profile become non-analytic for certain values of the input/output rates. This leads to a phase 
diagram with three main regions : a low density phase (when typically the input rates are small and the output rates 
are large), a high density phase (when the input rates are large and the output rates small), and a maximal current 
phase (in which the bulk density is 1/2 regardless of the boundary rates). 




RESERVOIR RESERVOIR 



FIG. 5: The Totally Asymmetric exclusion process on an open chain 

Although the ASEP with open boundaries was known to be an integrable system (Essler and Rittenberg 1996), 
the explicit diagonalisation of the Markov matrix by Bethe Ansatz could not be carried out because of the boundary 
terms that violate the conservation of the total number of particles. Only numerical and phenomenological studies of 
the spectrum were available (Bilstein and Wehefritz 1997, Dudziiiski and Schiitz 2000). In a remarkable recent paper, 
de Gier and Essler (2005) were able to derive the Bethe Ansatz equations describing the complete spectrum of the 
Markov matrix of the partially asymmetric exclusion process with open boundaries and general input/output rates. 
Their work uses an exact solution of the open XXZ chain with non-diagonal boundary terms (Nepomechie 2003 and 
2004, Cao et al. 2003). De Gier and Essler have calculated analytically the spectral gap of the TASEP with open 
boundaries and have discussed the various regimes throughout the phase diagram; in particular they have discovered 
a fine structure of the phase diagram that could not be found from the stationary state alone. Their work opens 
the way to many further investigations, such as the description of the excited states and the calculation of the large 
deviation function of the current in the open system. 

VII. CONCLUSION 

The Bethe Ansatz was initially used to diagonalize Quantum Hamiltonians (such as the Heisenberg Spin Chain, 
the one-dimensional Bose gas with 5-interactions) and to calculate real spectra of Hermitian Operators. In the 70's 
and the beginning of the 80's, the exact solutions of the six and eight vertex models and the development of the 
quantum inverse scattering method (Baxter 1982, Faddeev 1984) showed that many systems of interest in equilibrium 
statistical mechanics belonged to the class of integrable models. The ASEP, which plays the role of a paradigm in 
non-equilibrium statistical mechanics, thanks to its rich and complex phenomenology, has been shown to be integrable. 
This integrability property has allowed to derive many exact solutions for the ASEP and for some related models that 
we have reviewed here. These analytical results add to our knowledge of systems far from equilibrium and allow us 
to compare them with systems at equilibrium on a quantitative level. 

The use of integrability techniques in non-equilibrium statistical mechanics is an active field of current research. 
In particular, we believe that the following problems should lead to interesting results in near future : the use of the 
nested Bethe Ansatz for ASEP with multiple species of particles, the investigation of ASEP with open boundaries, the 
study of exclusion processes with disorder, and the application of integrability techniques to the zero-range process 
(Povolotsky 2004, Povolotsky and Mendes 2006) 
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